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PREFACE 
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reproduced  in  the  Appendix  for  technical  guidance  and  help  in  preparing  the  manu  x;ript  of 
the  paper. 
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INTRODUCTION 


Described  in  this  report  is  a  preliminary  model  for  predicting  the  performance  of  the 
MSX  SPIRIT  III  interferometer.  The  model  provides  estimates  of  the  expected  noise  and 
saturation  levels  of  the  instrument  in  measuring  spectral  radiances  of  the  Earth  and  its 
limb.  (The  Appendix  and  its  references  describe  the  the  MSX  program  and  SPIRIT  III 
interferometer). 

This  "noise  mcxiel"  of  the  SPIRIT  III  interferometer  is  implemented  as  a  computer 
code  module  named  IFR_NOIS,  which  will  be  part  of  the  Earthlimb  automated  analysis 
capability  for  processing  MSX  data  at  the  Phillips  Laboratory  Data  Analysis  Center  (PL 
DAC).  The  mcxJule's  primary  subroutine  is  IFR_NOISE.  In  the  DAC  automated  spectral 
analysis,  this  subroutine  (as  well  as  many  others  involved  in  screening,  indexing  and 
cataloging  the  interferometer  data)  will  be  "driven"  by  the  routine  IFR_SPC_DRIVER. 
The  current  version  of  IFR  NOISE  is  Revision  Ic,  dated  28  April  1994. 

The  primary  purpo.se  of  the  present  report  is  tti  provide  a  u.ser's  guide  for  subroutine 
lFR_NOISE  and  to  present  some  performance  predictions  ba.sed  on  this  noise  mixJel.  The 
first  three  sections  define  the  purpose  of  the  routine,  and  provide  related  background 
material  and  detailed  usage  iastruciions.  The.se  .sections  constitute  the  user's  guide. 

The  .section  entitled  Theoretical  Basis  fully  dcKuments  the  equations  and  conventions 
upon  which  the  noi.se  calculations  are  ba.scd.  it  a.s.sumes  familiarity  with  the  interfcn>meter 
system  design,  which  is  described  in  Ref  1  and  in  works  cited  in  Section  1.2  of  Ref.  1 .  A 
discu.ssion  of  the  qualitative  effects  of  the  .several  noise  sources  treated  in  the  present  noise 
mixJel  will  be  found  in  Ref  2.  (Also  given  in  Ref  2  are  interferometer  perfi>rmance 
predictions  ba.sed  on  an  earlier,  very  crude  noise  model;  these  are  superseded  by 
predictions  ha.sed  on  the  mixlel  de.scribed  herein.) 

The  section  entitled  Interferometer  Parameters  lists  the  iastpiment  specificatirms  u.seJ 
in  the  current  noi.se  mixlel,  such  as  the  fields  of  view,  collector  area,  etc.,  that  determine 
its  capabilities.  These  data  are  represented  as  constants  in  the  subroutine.  (Filter  spectral 
transmissions,  which  are  supplied  to  the  routine  via  common  storage,  are  not  listed  ) 
Some  t>f  the  "specification.s,"  such  as  the  interferometer's  wavelength-dependent 
mcxlulation  index,  are  ba.sed  on  engineering  data  and  recent  tests  at  SDL.  The  li.sted 
specification.s  are  believed  to  be  up-to-date  as  t>f  early  July,  199.1.  The  subroutine  is 
.structured  such  that  the  interferometer  parameter  values  are  easily  changed. 

The  remaining  .sections  give  the  channel  Noi.se  Equivalent  Spectral  Ry 
(NESRs)  and  other  predictions  of  the  m>ise  model.  The  Appendix  is  a  preprint  t>f 
paper  that  used  the  present  noi.se  mixlel  to  predict  the  performance  i>f  the  SPIK 
interferometer  in  planned  MSX  experiments  for  remote  sensing  of  atmospheric  trace 
constituents. 


It  is  emphasized  that  this  report  describes  a  preliminary  version  (Rev.  Ic)  of  the 
interferometer  noise  mcxJel,  The  final  version  of  the  nicxlel  (module  IFR_NOIS,  Rev.  2) 
will  be  delivered  to  PL  in  May  or  June  1994.  Revision  Ic  is  e.ssentially  a  predictive  model 
based  on  given  "efficiencies"  for  the  system's  optical  components  and  detectors.  Revision 
2  will  use  preflight  calibration  data  in  place  of  most  of  these  efficiercies;  it  will  be  more  of 
an  empirical  mtxlel.  The  present  report  .supersedes  the  ARC  technical  report  describing 
Revision  0  of  the  noise  model. ^ 


PURPOSE 

Subroutine  IFR_NOISL  can  be  u.sed  to  e.stimate  the  total  system  noi.se  for  the  .MSX 
SPIRIT  III  interferometer.  The  returned  total  noi.se  includes  the  effects  of 

a)  preamp-John.son  noi.se  (sy.stem  dark  noi.se,  which  determines  the  NESR), 

b)  photon  noise,  and 

c)  digitization  noise  from  A/D  conversion. 

The  spectral-domain,  total-noise  e.stimate  depends  on  the  length  and  apodization  of  the 
iiiterferogram  (i.e.,  on  the  effective  spectral  resolution)  and  on  the  .scene  spectral  radiance 
observed  by  the  interferometer.  The  .scene  contributes  photon  noi.se;  it  al.so  determines  the 
automatic  gain  setting  which,  in  turn,  determines  the  A/D  noi.se  level. 

The  subroutine  can  also  be  used  to  compute  just  the  Noi.se  Equivalent  Spectral 
Radiance  (NESR)  or  t('  return  just  a  scale  factor  SF  by  which  the  NESR  shi>uld  be 
multiplied  to  obtain  the  tc^tal  noise. 

BACKGROUND 

The  aultimated  analy.sis  of  MSX  SPIRIT  III  interferometer  data,  as  charted  in 
Fig.  .^-.^6  of  the  Earthlimb  Automated  Analysis  Plan  (EAAP)*^,  produces  a  template 
spectrum  and  noise  spectrum  for  each  measured  (CONVERT)  spectrum  that  is  analyzed. 
Summary  prtidiict  spectra  will  show'  the  measured  spectral  radiance  and  total  noise 
spectral  radiance  m  the  same  plot  (usually  the  larger  of  the  two  values  at  e;ich 
wavelength).  Subroutine  IFR  NOISE,  and  tne  section  of  the  calling  procedure  that 
assigns  the  actual  arguments  for  the  subroutine,  constitute  the  procedure  labeled 
"ESTIMATE  SYSTEM  TOTAL  NOISE  SPECTRUM"  in  the  EAAP  fiowchart. 

The  use  of  different  operational  mtxles  in  the  subroutine  offers  flexibility  in  how  noi.se 
estimates  are  obtained.  For  example,  the  calling  prcKedure  can  compute  NESRs  for  the 
six  channels  once-and-for-all.  aru!  thereafter  u.se  the  routine  only  to  obtain  scale  factors 
SE'.  It  can  elect  to  use  onh  NESRs  obtained  from  engineering  test  data  or  a  DCATT 
calibration  sequence  and  scale  ihe.se  by  SF  to  obtain  the  total  noise. 


The  argurrients  supplied  by  the  calling  procedure  include  interferogram 
iength/apodization  parameters  and  a  "source"  spectrum.  The  source  spectrum  (which  is 
required  only  when  total  noise  or  SF  is  to  be  computed)  will  usually  be  one  obtained  from 
measured  data  via  the  Canonical  CONVERT  program;  the  Iength/apodization  parameters 
would  be  those  u.sed  in  the  CONVERT  processing.  For  modeling  expected  performance 
or  when  the  S/N  in  the  observed  spectmm  is  expected  to  be  low,  the  calling  procedure 
may  want  to  supply  a  template  spectrum  as  the  source  .spectrum.  The  calling  procedure 
can  either  supply  a  gain  setting  or  direct  the  routine  to  predict  the  automatic  gain  level. 


USAGE 

Subroutine  IFR_NOISE  requires  spectral  tran.smittance  data  for  the  six  interferometer 
filters;  this  data  is  assumed  to  exi.st  in  the  two  dimensional  array  filtjrans  in  common 
block  GETJFR  DAT.  The  subroutine  is  invoked  by  the  Fortran  statement 

CALL  IFR_NOISE(  MODE,  ICHAN,  LASER,  lATYP,  ALPHA,  LEN_APO, 

RESLN,  GAINACT,  GAIN2,  GAINl,  SPECTRUM,  LEN_SPC,  WAVNl,  DELW, 
TOT_NOIS,  SF.  BPR,  ICLIP,  !ERR.  SATBPR,  EFFM  ) 

where: 

MODE  determines  the  operational  m«xJe; 

=  - 1  if  the  NESR  is  to  be  returned. 

=  0  if  the  scale  factor  SF  ({total  noi.se  }/NESR)  is  to  be  returned. 

=  1  if  both  total  noise  (NESR  times  SF)  and  SF  are  to  be  returned. 

The  NESR  is  returned  in  array  SPECTRUM.  Other  parameters  (see  below)  determine 
its  spectral  range  and  the  spacing  of  the  spectrum  samples. 

ICFMN  =  the  channel  number  for  the  NESR,  SF,  or  total  noise. 

LASER  determines  which  laser  reference  is  in  operation: 

=  (t  if  the  primary  (HeNe)  laser  is  in  use 
=  1  if  the  solid  slate  l.iser  is  in  use 

lATYP  selects  the  apodizalion  function  to  be  represented  in  finding  the 
spectral  resolutiim  that  the  NESR  or  total  noise  will  correspond  to: 

=  0  for  no  apodizalion 

=  I  for  triangular  apodizalion 

=  2  for  Kaiser-Bes.sel  (K-B)  apodizalion 

Al.PHA  =:  value  of  parameter  n  in  the  K-B  window  function  (if  lATYP  -  2) 

l^E'N  APO  =  number  vif  samples,  after  two-sided  Irunc.ttion 'apodizalion  to  be  used  in 
finding  the  spectral  resolution. 
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RESLN  =  the  computed  spectral  resolution  (returned  in  cm'') 


GAINACT,  if  not  zero,  is  a  supplied  gain  for  channel  ICHAN. 

If  the  supplied  GAINACT  is  zero,  the  routine  sets  the  gain  internally 

to  the  computed  value  GAIN2.  (When  MODE  =  -1,  the  supplied  GAINACT 

IS  not  used  and  GAIN2  is  not  computed/returned.) 

GAIN2  =  the  theoretical  automatic  gain  as  determined  from  the  supplied  spectral 
radiance  SPECTRUM,  and  assuming  the  band  photon  radiance  is  increasing 
(high-gain  side  of  hysteresis  loop);  GAIN2  is  computed  and  returned  if 
MODE  is  not  -1  (whether  or  not  GAINACT  is  zero). 

GAINl  has  the  same  definition  as  GAIN2,  except  decreasing  band  photon  radiance 
is  assumed  (low-gain  side  of  hysteresis  loop). 

SPECTRUM  is  a  supplied  or  returned  array  of  spectral  radiances: 

If  MODE  is  not  -1,  it  is  a  supplied  radiance  spectrum  over  channel  ICHAN 
representing  what  the  interferometer  "sees,"  i.e.,  a  spectrum  from 
Canonical  CONVERT  or  a  predicted  (template)  spectrum.  Parameters 
LEN_SPC,  WAVNl  and  DELW  define  the  spectral  range  and  spacing.  Values 
in  a  supplied  SPECTRUM  are  expected  in  the  units  W/cm^  sr  cm"'. 

If  MODE  is  -1,  it  is  the  returned  NESR  spectrum  in  W/cm^  sr  cm"'. 

LEN_SPC  =  the  length  of  array  SPECTRUM;  LEN_SPC,  WAVNl  and  DELW  are 
supplied  (input)  parameters. 

WAVNl  =  the  wavenumber  (cm"')  represented  by  the  first  element  of  SPECTRUM. 

DELW  =  the  wavenumber  spacing  (cm*')  of  spectral  radiances  in  SPECTRUM. 

1  OT_NOIS  =  the  returned  total  noise  spectrum  in  W/cm^  srcin''  if  MODE  =  1. 
LEN_SPC,  WAVNl  and  DELW  define  the  spectral  range  and  spacing.  When 
MODE  is  -1  or  0,  array  TOT_NOIS  is  u.sed  as  workspace;  i.e.,  its  contents  prior 
to  the  .subroutine  call  are  not  preserved. 

SF  (if  MODE  is  0  or  1)  is  the  returned  .scale  factor.  Multiplying  a  NESR 
spectrum  by  SF  converts  it  to  a  total  noise  spectrum. 

BPR  =  the  Band  Photon  Radiance  (ph/s  cm^  sr)  corresponding  to  SPECTRUM 
if  MODE  is  0  or  1  (.see  Theoretical  Basis  Section). 

ICLIP  is  returned,  if  MODE  is  Oor  1,  to  indicate  whether  the  interferometer 
channel  is  predicted  to  be  clipped; 

=  0  if  the  fringes  are  not  clipped  fo:  the  supplied  SPECTRUM 
=  !  if  the  fringes  are  dipped 
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lERR  is  a  returned  error  indicator: 

=  0  if  no  error 

=  1  if  the  spectral  filter  for  channel  ICHAN  is  not  completely  within  the  spectral 
range  defined  by  the  supplied  values  of  LEN_SPC,  WAVNl  and  DELW. 
(Filter  bands  are  defined  by  arrays  lamls  and  lam2s\  this  data  is  assumed  to 
exist  in  common  block  FILT  CUTOFFS.) 

SATBPR  -  the  BPR  that  causes  interferogram  central  fringe  clipping;  returned  only  if 
MODE  is  0  or  1  (see  BPR  and  Theoretical  Basis  Section). 

EFFM  =  the  effective  modulation  index,  defined  as  the  computed  central  fringe  ac 

interferogram  value  divided  by  the  computed  dc  level;  returned  only  if  MODE  is 
0  or  1 . 

Note  that  the  calling  procedure,  which  is  normally  part  of  the  automated  SPIRIT  III 
spectral  analysis,  will  have  knowledge  of  how  individual  interferograms  were  processed  by 
Canonical  CONVERT.  In  other  words,  the  calling  procedure  is  supplied  values  for 
ICHAN,  LASER,  lATYP,  ALPHA  and  LEN_APO  by  Standard  CONVERT  and/or 
Canonical  CONVERT  and  the  DCE  Experiment  Script.  The  calling  procedure,  in  turn, 
supplies  these  values  to  subroutine  IFR_NOISE. 

The  largest  values  to  be  expected  for  LEN_APO  are  roughly  8220  and  4110  samples, 
respectively,  for  even  and  odd  ICHAN  values.  These  correspond  to  the  maximum  scan 
time  of  4.3  s,  with  0.14  s  lost  to  turnaround.  The  total  retardation  change  X  over  the  two- 
sided  .scan  equals  drive  velocity  times  effective  scan  time,  approximately  0.25  OPD  cm/s 
times  4.16  s,  or  1.04  cm  (OPD  is  optical  path  difference,  or  retardation).  The  above 
estimates  for  LEN_APO  equal  X/Ax,  where  Ax  is  one  of  the  two  HeNe  sampling  intervals 
given  by  Eq.(l)  in  the  Theoretical  Basis  section.  Use  of  these  values  for  LEN_APO 
implies  that  the  long  interferograms  have  been  zero-filled  to  the  transform  array  sizes 
16384  and  8192,  respectively. 

An  anticipated  CONVERT  processing  option  is  to  slightly  truncate  the  long 
interferograms  to  8192  and  4096  .samples,  in  which  case  these  would  be  the  appropriate 
values  to  supply  for  LEN  APO.  The  corresponding  unapixlized  resolution  is 
approximately  0.96  cm'l;  this  is  the  value  of  RESLN  that  would  be  returned  when 
lATYP  =  0. 


THEORETICAL  BASIS 

This  section  defines  the  noise  mcxlel  for  the  SPIRIT  III  interferometer  in  terms  of 
equations  and  conventions  implemented  in  subroutine  IFR_NOISE.  The  subroutine's 
dummy  arguments  appear  in  capital  letters  and  are  as  defined  in  the  above  Usage  section. 
(Identical  nt>tation  is  used  in  this  section,  in  the  above  Usage  section  and  in  the  Fortran 
source  code  for  the  dummy  arguments,  but  not  for  other  quantities.)  The  order  in  which 
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quantities  are  calculated  in  the  subroutine  and  in  this  description  of  the  noise  model  is 
roughly  the  same. 

The  input  parameter  LASER  is  used  to  specify  which  laser  is  being  used  in  the 
interferometer  reference  channel  to  sample  the  interferogram.  If  LASER  is  zero,  the 
primary  (HeNe)  laser  is  in  use,  and  the  laser  wavelength  is  assigned  the  value 
0.632985  pm.  If  LASER  is  not  zero,  the  .solid  state  laser  is  being  u.sed,  and  set  to 
0.6700  pm. 

The  input  parameter  ICHAN  is  the  interferometer  channel  number  for  the  current 
subroutine  calculations.  The  .sampling  interval  Ax  is  twice  Xl  for  the  even-numbered 
channels  and  four  times  Xl  for  the  odd-numbered  channels;  i.e.,  the  even-numbered 
channels  are  located  at  relatively  short  wavelengths  (high  wavenumbers)  and  mu.st  be 
sampled  at  a  higher  rate  than  the  odd-numbered  channels.  The  subroutine  obtains  Ax  in 
cm  as 

Ax  =  (MOD(ICHAN,2)+ 1  ).2.0*XL*10-f^.  (1) 

The  calculated  Ax  and  other  quantities  are  used  to  estimate  the  spectral  resolution  and  the 
parameter  kl,  defined  below. 

The  length,  in  .samples,  of  the  two-sided  interferogram  after  truncation,  if  any 
(including  any  truncation  effected  by  an  apodization  that  the  calling  procedure  wants  to 
represent),  is  specified  by  the  input  parameter  LEN_APO.  The  corresponding  length  in 
cm  is  Ax  times  LEN_APO.  The  reciprocal  of  this  length  is  the  resolved  spectral  interval 
Av,  if  the  interferogram  was  not  apodized.  The  subroutine  calculates  a  factor  rf{resfac  in 
the  source  code)  that  accounts  for  the  effect  on  re.solution  of  any  apodization  that  the 
calling  procedure  wants  to  be  repre.sented  in  the  noi.se  calculation.  Thus,  the  re.solved 
spectral  interval  in  cm'^  is  obtained  as 

Av  =  r//(  Ax.LEN_APO  )  (2) 

The  subroutine  allows  the  specification  of  three  apodization  "types"  via  the  input 
parameter  lATYP;  Values  of  0,  1  and  2  corre.spond  to  "no  apodization,"  triangular 
apodization  and  Kaiser-Bessel  (K-B)  apodization,  respectively;  the  length  of  the 
apodization  window  is  LEN_APO.  When  lATYP  is  2,  the  additional  input  parameter 
ALPHA  .specifies  the  value  of  a  in  the  K-B  window  function,  where  the  vv'indow  function 
and  a  are  as  defined  by  F.  J.  Harris.^  The  factor  rf  evaluates  to  1 .0  if  lATYP  =  0,  to  2.0  if 
lATYP  =  1,  and  to  ALPHA  if  lATYP  =  2.  The  computed  Av  is  the  distance  to  the  first 
zero  in  the  instrument  line  shape  (which  is  the  Fourier  cosine  transform  of  the  apodization 
window  function).  The  computed  Av  is  exact  when  lATYP  is  0  or  1,  and  is  approximate 
when  lATYP  is  2.  The  subroutine  returns  Av  as  RESLN. 
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7'he  parameter  kl  calculated  by  the  routine  converts  quantities  normally  specified  in 
per  root-Hz  units  (such  as  NEP)  to  per  wavenumber  (cm'*)  units.  As  shown  by  Wyatt^, 
this  parameter  can  be  cast  in  the  form 

kl  =  (vo/Av)l/2,  (3) 

where  vq  is  the  drive  velocity  of  the  interferometer  in  OPD  cm/s.  The  calculated  kl  has 
the  units  cm  jhe  calculated  NESR  and  total  noise  includes  kl  as  a  factor;  i.e., 

these  noise  spectra  vary  as  Av  '*''2. 

The  calculation  of  the  interferometer  noise  requires  the  spectral  transmission  of  the 
filter  in  channel  ICHAN.  This  spectrum  is  supplied  to  the  routine  via  an  array  in  common. 
Hereafter,  the  filter  spectral  transmission  for  channel  ICHAN  will  be  denoted  tp(i),  where  i 
is  a  spectrum  index.  For  channel  5,  tp(i)  includes  the  effect  of  IFOV  reduction  by  a  mask. 

In  preparation  for  the  calculation  of  spectrum  noise,  the  routine  obtains  the  mean- 
square  preamp  noi.se  voltage  Vp,  the  mean-square  thermal  noise  voltage  vj  and  the 
resultant  rms  dark  noise  voltage  Vjjf^: 


Vp  =  (27iRLCTEn)2F3/3, 

(4) 

V,  =  4kTRLF, 

(5) 

vdn  =  (Vp  +  v,  )*/2, 

(6) 

where 


Rl  is  the  feedback  resistance  of  the  Traas-Impedance  Amplifier  (TIA)  preamp, 

Cj  is  the  TIA  input  capacitance, 

Ep  is  the  TIA  input  noi.se  voltage, 

T  is  the  temperature  of  the  preamp  and  detectors,  and 
F  is  the  information  noise  bandwidth  of  channel  ICHAN. 

The.se,  of  course,  are  voltages  in  the  time  (interferogram)  domain. 

The  noise-equivalent  spectral  radiance  (NESR)  depends  on  the  above  dark  noise 
voltage  V(jn,  various  constant  efficiency  factors  and  efficiencies  dependent  on  channel 
and/or  wavelength,  like  beamsplitter  net  efficiency,  mcxJulation  index,  tp(i)  and  detector 
quantum  efficiency.  The  routine  has  a  kxip  over  wavenumber  that  searches  and 
interpolates  in  pre-Jefined  data  arrays  to  establish  values  for  most  of  the  spectrally  varying 
quantities.  In  the  same  loop  the  routine  integrates  over  an  optionally  supplied  scene 
radiance  spectrum  to  obtain  the  "Band  Photon  Radiance"  (BPR).  The  BPR  determines  the 
photon  noise  contribution  to  the  total  noi.se,  as  well  as  the  mtiximum  of  the  ac  component 
of  the  interferogram.  The  latter  is  used  to  estimate  the  automatic  gain  setting  and  the 
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resultant  A/D  noise  contribution.  The  spectral  integration  art!  the  calculation  of  BPR  and 
gain  is  not  performed  when  the  calling  procedure  specifies  that  only  the  NESR  is  to  be 
returned, 

The  calling  procedure  specifies  a  wavenumber  mesh  and  spectral  range  over  which  the 
NESK  or  total  noise  will  be  calculated.  This  is  also  the  mesh  on  which  the  scene  spectral 
radiance  is  supplied  if  the  total  noise  or  .scale  factor  SF  is  to  be  calculated.  The  routine 
finds  the  detector  quantum  efficiency  q(i)  (as.sumed  to  be  the  same  for  all 
detectors/channehs)  and  the  filter  spectral  transmission  tp(i)  interpolated  onto  this  mesh. 
The  vectors  q(i)  and  tp(i)  will  now  denote  interpolated  values  rather  than  values  in  a  data 
statement  or  common  array;  these  are  indexed  by  wavenumber,  whereas  the  original  data 
values  are  defined  on  a  wavelength  grid. 

Let  X(i)  denote  the  wavelength  in  p.m  corresponding  to  index  i,  i.e.,  to  wavenumber 
v(i)  m  the  user-defined  mesh.  The  '  oltage  responsivity  r(i)  of  the  detector  at  this 
wavelength  is  calculated  as 

KO  =  r20*(  M')/20  )*q(i)/q20.  O) 

where  r2o  and  q2o  are  the  detector  responsivity  and  quantum  efficiency,  respectively,  at 
wavelength  20  pm.  A  corresponding  NEP  is  obtained  as 

NEP(i)  =  vdn/(Fl/2*r(i)).  (8) 

Finally,  the  NESR(i)  is  calculated  as 

NESR(i)  =  kl*NEP(i)/(  Ac*Q*te*tp(i)).  (9) 

where 

Ac  is  the  collecting  area  of  the  entrance  aperture, 

Q  is  the  .solid  angle  IFOV  for  channel  ICHAN, 

^e  ~  U*lbs*’^ 

t(  is  the  telescope  reflectivity  for  channel  ICHAN, 

tj,jj  is  the  FTS  cube  (beamsplitter,  etc.)  net  efficiency,  and 

m  is  the  modulation  index. 

Both  tj^s  and  m  are  functions  of  wavelength  in  the  current  version  of  the  noise  model.  The 
NESR  values  returned  by  the  subroutine  have  the  units  Watt.s/cm2  sr  cm'*. 

An  alternative  form  for  Eq.  (7)  is 


r(i)  =  r2o*Rd('). 


(10) 


where 


Rd(i)  =  (X(i)/20)*q(i)/q2o  (11) 

is  the  detector  relative  response.  Similarly,  Eq.  (8)  can  be  wriiten  as 

NEP(i)  =  NEP(Xp)  /  Rd(i),  (12) 

where  NEP(>Lp)  =  vdn/(  Ih®  NEP  at  p>eak  spectral  response,  i.e.,  the 

response  at  X  =  Xp  «  20  pm.  Substituting  Eq.  (12)  into  (9)  and  writing  Rd(i)*tp(i)  =  Rs(i), 
where  Rs(i)  is  the  "system  spectra!  response  function,"  results  in  the  following  alternative 
expression  for  the  NESR: 

NESR(i)  =  kl*NEP(Xp)/(  Ac*Q*te*Rs(i)).  (13) 

USU  provided  a  file  of  the  detector  relative  response  Rj.  This  file  and  Eq.(ll)  were 
used  to  obtain  a  set  of  quantum  efficiencies  q  for  subroutine  IFR_NOISE. 

As  indicated  earlier,  a  user-supplied  radiance  spectrum  is  required  only  when  the  total 
noi.se  and/or  scale  factor  3F  is  to  be  returned.  The  supplied  spectrum  will  be  denoted 
Ny.  It  represents  spectral  radiance  in  per  cm’*  units  versus  wavenumber  v.  The 
corresponding  interferogram  voltage  vs.  retardation  x  at  the  TIA  is 


vs(x)  =  R,  e  t  A  Q 
^  ^  L  t  c 


i 


N  q  t„  t. 
v^v  Fv  bs 


[1  +  mcos(2n;vx)] 


hcv 


dv, 


Chan 


(14) 


where  e  is  the  elementary  charge  and  he  is  Planck's  constant  times  the  speed  of  light.  In 
this  equation,  qy  and  tpy  represent  the  quantities  q(i)  and  tp(i),  defined  earlier.  The 
integral  is  over  the  wavenumber  range  of  channel  ICHAN. 

The  Band  Photon  Radiance  (BPR)  is  defined  as  the  effective  dc  photon  radiance 
(ph/s  cm2  sr)  at  the  entrance  aperture: 

C  N  q  tj,  t, 

BPR=  [  Y.Fy,bs 

J  hcv 

Chan 

Analogous  to  the  dc  and  ac  components  of  the  interferogram  (see  Eq.  14),  we  define  BPR 
to  be  a  dc  term  and  BPRm  to  be  the  corresponding  ac  term: 
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1*  N  q  t  rr  •  L. 

BPRm  =  I  m  -  ^  dv.  (15b) 

J  hcv 

Chan 

Hereafter,  vs3^;  will  denote  the  ac  voltage  at  the  central  fringe  (at  x  -  0)  and  vsjj.  will 


denote  the  dc  voltage.  Using  Eqs.(14)  and  (15)  these  voltages  can  be  written 

v.S(jc  =  RL*c*(t|/2)»A^»Q*BPR,  (16a) 

''■‘’ac  =  RL*e*(tt/2)*Ac*0»BPRm.  (16b) 

The  dc  photon  spectrum  (photons/s  cm'*)  impinging  on  the  detector  is 

<tpv  =  (tt/2)*A<;.*Q*Ny*tpy*t;^/hcv.  (17) 

Using  Eqs.(15a)  and  (17)  we  can  write  (16a)  in  the  form 

”dc  =  ■'l'  J'l’pv'iv 

Chan 

The  rms  of  the  photon  noise  voltage  throughout  the  iriterferogram  is  given  by 

"pn  =  (19) 

Chan 

which,  from  Eq.(18),  can  be  written 

vpn  =  (RL*e*vsdc*2F)l/2.  (20) 


The  routine  obtains  BPR  and  BPRm  according  to  Eqs.(15),  vsgc  and  vsj^.  from  Eqs.(16), 
and  then  calculates  Vpp  from  Eq.(20).  Yet  to  be  determined  is  the  A/D  noise  level  in  the 
interferogram. 

The  number  of  counts  C  out  of  the  12-bit  AT)  converter  corresponding  to  a  given 
interferogram  ac  maximum  voltage  vsg^  is 
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('  =  I  (  - : ) .  (  vv,^.  4  5  V  )  1(1  V  )  4  1  ) .  Cl. 


Cl) 


whfft*  Cl  is  thf  gum  selling  Thus,  iht  number  i>f  counis  prtHJu*.ed  al  uiiil  gam  by  ihe 
svsiem  when  vS;,,^^  is  -5.  0  and  5  vtWis  is  1.  204^^  and  4(1*^.''  (=  2*‘  •  1).  respectively  In 
aulomalic-gain  miKle.  the  gain  is  decrea.sed  by  a  factor  of  two  each  lime  an  incri’ustnfi 
number  of  counts  a:  the  interferogram  m;wimum  exceeds  a  ceiiain  C'  value.  Similarly,  a 
faclor-of-two  gain  increase  ikcuis  when  the  maximum  counts  are  di\rca.\tnf;  and  falls 
beh>w  a  particular  C  value.  At  unit  gain  the  two  voltages  that  effect  auti'matic  gain 
changes  are  3.7.*'  v.  and  3.7,*'/4  =  0.^37.*'  v.  For  an  arbitrary  gain  setting  G  we  conclude 
(using  the.se  two  voltage  levels  and  Eq  21)  that  gain  changes  occur  when 
|C-2048|  equals  2*’’  or  2*’*,  where  bl  =  8.58426  and  f2  =  10.58426.  This  is  the  same  as 
saying  automatic  gain  changes  correspond  to  an  irierferogram  ac  maximum  (absolute 
value)  which  is  equivalent  to  approximately  8.584  or  10.584  bits  at  the  current  gain 
setting. 

Subroutine  IFR_NOISE  uses  Eq.(21),  the  above  bl  and  b2  values,  and  the  known 
maximum  gains  for  the  six  channels  to  determine  the  automatic  gain  setting  corresponding 
lo  vsg^..  The  routine  actually  returns  two  values,  GA1N2  and  GAINl,  corresponding  to  an 
increasing  and  decreasing  signal,  rj.spectively.  As  de.scribed  below,  it  al.so  determines 
whether  the  interferogram  central  fringes  are  predicted  to  be  clipped. 

The  calling  procedure  can  supply  a  value  for  input  parameter  GAINACT,  which  if 
non-zero  is  used  by  the  routine  as  the  actual  gain  .setting  (which  has  been  reported, 
presumably,  by  CONVERT).  If  the  .supplied  value  is  zero,  the  routine  uses  the  calculated 
GAIN2  as  the  current  gain  setting.  The  maximum  gains  Gn^g^  for  the  six  interferometer 
channels  were  supplied  by  USU,  and  are  defined  'n  a  DATA  statement  in  the  subroutine. 

The  A/D  noi.se  voltage  at  the  TIA,  which  will  be  denoted  v-jj,  is  defined  as  that 
corresponding  to  an  uncertainty  of  1/2  in  the  count  level  C.  By  using  Eq.(21)  to  equate  a 
change  AC  of  1/2  to  a  change  Avs^^;  =  v^j  and  .solving  for  V3J,  we  obtain 

vad  »  (  l()v./2l2)/(2.G),  (22) 

where  G  is  tlie  current  gain  setting.  This  equation  is  used  by  the  routine  to  obtain  the  A/D 
noi.se  voltage. 

The  routine  obtains  the  total  noise  voltage  at  the  TIA  as  the  root -sum -square  of  the 
three  noise  components: 


tn 


=  ( 


dn 


+  V 


pn 


(23) 


Finally,  it  obtains  the  returned  scale  factor  SF  as  the  ratio  of  total  noise  to  dark  noise: 


SF  =  v,p  /vjp. 


(24) 


II 


Although  SF  is  obtained  from  voltages  in  the  time  domain,  it  is  also  the  ratio  of  total 
spectrum  noise  to  the  NESR.  Thus,  the  total  noise,  which  is  returned  as  the  array 
TOT  NOIS,  is  calculated  fn>m 

TOT_NOIS(i)  ^  SF*NESR(i).  (25) 

The  gain  can  have  one  of  seven  levels  ranging  between  and  =  G^j,^/2^. 
Clipping  of  the  interferogram  occurs  when  the  gain  has  its  minimum  value  while  |C-2048| 
exceeds  10* The  routine  finds  the  correspcmding  BPR  level  above  which  the 
interferogram  is  clipped  and  returns  it  as  SATBPR.  The  routine  also  returns  a  value  for 
ICLIP;  it  is  I  if  BPR  i  SATBPR  and  0  otheswise. 

The  routine  also  calculates  and  returns  an  "effective"  modulation  efficiency  EFFM, 
defined  as  the  ratio  of  ac  to  dc  voltage  at  the  interferogram  center,  i.e.,  EFFM  =  vSacA'Sj^, 
=  BPRm  /  BPR. 

Many  interferometer  parameters,  such  as  A^  and  the  wavelength-dependent 
beamsplitter  efficiency  tj^,  are  assigned  values  through  PARAMETER  and  DATA 
statements  in  the  subroutine  code.  This  is  also  true  of  the  full-scale  A/D  i.nput  voltage  V{j;, 
the  number  of  A/D  bits  Ngj,  the  gain-switching  bit  levels,  etc.,  even  though  their  assigned 
(actuc!  SPIRIT  III)  values  have  been  used  in  the  above  text  and  equations  to  simplify  the 
discussion;  e.g.,  vj^  =  10  v.  and  Ngj  =  12. 

The  wavelength-dependent  mcxlulation  index  m  is  computed  from  a  fifth-order 
polynomial  evaluated  by  a  short  function  routine  included  in  module  IFR_NOIS. 
Subroutine  IFR_NOISE  includes  a  data  statement  defining  the  six  coefficients  needed  by 
the  function.  Each  reference  to  the  function  routine  supplies  these  coefficients  and  defines 
the  wavelength  for  which  the  mcxlulation  index  is  to  be  computed. 


INTERFEROMETER  PARAMETERS 

This  section  lists  the  constants  in  subroutine  IFR_NOISE  that  define  the 
interferometer  noise.  Table  !  shows  only  those  con.stants  needed  by  the  routine  in 
calculating  the  NESR  corresponding  to  a  particular  spectral  resolution.  The  mcxlulation 
index  m  at  wavelength  W  is  computed  by  the  routine  as 

m  =  (  cl  +  W  (c2  +  W  (c3  +  W  (c4  +  W  (c5  +  W  c6) ))))  / 100,  with 
cl  =  2.734896,  c2  = -0.5068297,  c3  =  0.9076166,  c4  = -0.07671653 

c5  =  .002503216,  c6  =  -  0.00002928215.  (26) 


12 


Table  1.  Parameters  of  the  SPIRIT  III  interferometer  noise  mcxlel  that  affect  NESR 


Large  detector  IFOV  (sr) . 

Small  detector  IFOV  (sr) . 

Collector  area  (cm‘) . 

Drive  velocity  (OPD  cm/s) . 

PARAHETER/CHANNEL  1 

2 

3 

4 

.  1.37E-05 

.  5.14E-06 

.  393.0 

.  0.250 

5  6 

Telescope  reflectance...  0.932 
Information  noise 

0.911 

0.932 

0.922 

0.919  0.922 

bandwidth  (Hz) .  300 

1200 

600 

1200 

600  1200 

Filter  spectr'l  transT^ttancea :  per  files  supplied  on  floppy  disk 

by  USU/SDL  as  of  C/21/93 


TIA  feedback  resistance  (ohms) .  9.0E-f07 

TIA  input  capacitance  (farads) .  4.0E-11 

TIA  input  noise  (volts) .  l.OE-07 

TIA  operating  temperature  (deg  K) .  10 

Detector  responaivity  at  20  |im  (V/W) . 1.26E+09 


Detector  quantum  efficiency  Q  vs.  wavelength  k  in  )un: 


k 

0 

k 

Q 

k 

Q 

k 

Q 

2.0 

.4230 

2.5 

.3800 

3.0 

.2800 

3.5 

.2800 

4.0 

.3700 

4.5 

.4000 

5.0 

.3800 

5.5 

.3100 

6.0 

.3500 

6.5 

.3420 

7.0 

.3320 

7.5 

.3520 

8.0 

.3970 

8.5 

.3860 

9.0 

.4470 

9.5 

.4450 

10.0 

.5080 

10.5 

.5340 

11.0 

.5770 

11.5 

.6350 

12.0 

.6170 

12.5 

.6180 

13.0 

.6830 

13.5 

.7050 

14.0 

.7400 

14.5 

.6780 

15.0 

.7900 

15.5 

.7370 

16.0 

.9250 

16.5 

1.0300 

17.0 

.9330 

17.5 

.9060 

18.0 

.8460 

18.5 

.8120 

19.0 

.7900 

19.5 

.7810 

20.0 

.7930 

20.5 

.7740 

21.0 

.7960 

21.5 

.8660 

22.0 

.8650 

22.5 

.8650 

23.0 

.8650 

23.5 

.9000 

24.0 

.8640 

24.5 

.7770 

25.0 

.6600 

25.5 

.5640 

26.0 

.4880 

26.5 

.3670 

27.0 

.2660 

27.5 

.  1540 

28.0 

.0830 

28.5 

.0370 

29.0 

.0146 

29.5 

.0108 

30.0 

.0070 

30.5 

.0000 

Beeunsplitter  net 

efficiency  tbs  vs. 

wavelength  k  in 

^m: 

k 

tbs 

k 

tbs 

k 

tbs 

k 

tbs 

2.0 

.0250 

2.5 

.0265 

3.0 

.0280 

3.5 

.0420 

4.0 

.0600 

4.5 

.0650 

5.0 

.0700 

5.5 

.0730 

6.0 

.0760 

6.5 

.0840 

7.0 

.0960 

7.5 

.  1020 

8.0 

.1160 

8.5 

.  1100 

9.0 

.1050 

9.5 

.1160 

10.0 

.1270 

10.5 

.  1160 

11.0 

.1010 

11.5 

.  1000 

12.0 

.0980 

12.5 

.1200 

13.0 

.1400 

13.5 

.  1500 

14.0 

.1600 

14.5 

.1600 

15.0 

.1600 

15.5 

.  1600 

16.0 

.1670 

16.5 

.  1610 

17.0 

.1560 

17.5 

.  1560 

18.0 

.1560 

18.5 

.1470 

19.0 

.1370 

19.5 

.  1310 

20.0 

.1240 

20.5 

.  1180 

21.0 

.1120 

21.5 

.  1000 

22.0 

.0880 

22.5 

.0790 

23.0 

.0700 

23.5 

.0650 

24.0 

.0600 

24.5 

.0500 

25.0 

.0400 

25.5 

.0310 

26.0 

.0210 

26.5 

.0160 

27.0 

.0100 

27.5 

.0080 

28.0 

.0060 

28.5 

.0030 

29.0 

.0010 

29.5 

.0010 

30.0 

.0001 

30.5 

.0000 
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For  the  listeo'  value  of  drive  velocity,  the  internally  computed  frequency  conversion  factor 
kl  (see  Eq.  3  in  the  Theoretical  Basis  section)  evaluates  to  0.509  cm  for  the  case  of 
an  unaptxlized,  long  interferogram.  The  spectrally  resolved  interval  for  this  case  (see 
Eq.  2)  is  approximately  0.96  cm  *  if  the  interferogram  is  not  aptxiized. 

Table  2  lists  the  additional  constants  needed  to  predict  the  gain  level  and/or  total  noise 
for  a  given  input  radiance  spectnim. 

Table  2.  Additional  parameters  of  the  SPIRIT  III  interferometer  noise  model 


A/D  full-scale  input  volts .  10 

Number  of  A/O  converter  bits .  12 

AC  bits  for  ^ain  reduction . 10.58426 

AC  bits  for  gain  increase.. . 8.58426 

Maximum  number  of  gain  changes  allowed .  6 

Factor  by  which  gain  is  changed .  2.0 

PARAMETER/CHANNEL  12  3456 


Maximum  gain .  357.3  64.0  140. Q  64.0  140.8  64.0 


CHANNEL  NESRs 

A  simple  driver  program  was  written  to  test  subroutine  IFR_NOISE  for  MODE  =  -1 
and  obtain  its  prediction  of  the  six  channel  NESRs.  The  resultant  NESRs,  corresponding 
to  a  hmg,  unaptxlized  interferogram  (Av  =  0.96  cm"')  are  shown  in  Figure  1.  For 


Figure  1.  NESRs  predicted  for  the  SPIRIT  III  interferometer  by  the  noi.se  model 
code.  The  NESRs  correspond  to  a  spectral  re.solution  of  0.96  cm"',  unapodized. 
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comparison,  the  NESRs  measured  by  SDL  during  a  cold  lest  in  1993  are  shown  in 
Figure  2.  It  is  seen  that  the  noise  mtxlel  does  reasonably  well  in  reprtxlucing  the 
spectral  shapes  of  the  measured  NESRs,  but  does  not  accurately  predict  the  actual 
NESR  levels  for  all  channels.  The  mcxlel  coasi.stenily  undere.stimates  the  NESR  and 
predicts  substantially  lower  NESRs  than  are  measured  for  channels  2,  3  and  4. 


Figure  2.  NESRs  measured  by  SDL  (after  extrapolation  to  .spectral  resolution 
0.96  cm'*  and  removal  of  the  A/D  noi.se  component). 


These  discrepancies  led  to  the  decision  to  develop  another  revision  of  the  noise 
mcxlel.  As  stated  earlier,  the  forthcoming  Rev.  2  will  be  more  of  an  empirical 
mcxlel,  which  will  depend  on  measured  rather  than  predicted  NESRs.  It  will  still  be 
nec:e.ssary  in  the  Rev.  2  mcxlel  to  use  e.stimated  efficiencies  for  some  optical 
components  in  order  to  predict  the  A/D  and  photon  noise  contributions  to  the  total 
noise  spectrum.  However,  it  is  anticipated  that  the  revised  noise  mcxlel  will  provide 
much  better  accuracy  in  predicting  total  noise  levels  than  the  current  Rev.  Ic  model. 

The  current  interferometer  noi.se  mcxlel  was  developed  both  as  an  element  of  the 
MSX  automated  data  prcx;essing  and  as  a  tcx>l  for  estimating  the  performance  to  be 
expected  in  planned  MSX  experiments.  The  performance  predictions  described  in 
the  following  section  and  in  the  Appendix  are  ba.sed  on  a  temporary  mcxlification  of 
the  Rev.  Ic  mcxlel:  Predicted  NESRs  for  the  six  channels  are  each  scaled  by  a 
con.stant  .so  as  to  be  in  approximate  agreement  with  the  measured  NESRs  shown  in 
Figure  2.  The  scaled  NESRs  arc  then  combined  with  the  mcxlel's  prediction  of  A/D 
and  photon  noi.se  to  obtain  the  total  noise  spectrum.  The  Rev.  Ic  mcxlel  is  a.ssumed 
to  correctly  predict  gain  settings  and  the  cx:currence  of  clipped  interferograms. 


15 


OTHER  MODEL  PREDICTIONS 


Another  driver  ccxle  for  subroutine  GET  NOISE  was  written  to  obtain  estimates  of 
the  tangent  height  at  which  the  interferogram  becomes  clipped,  and  the  one  at  which  the 
maximum  S/N  in  the  spectrum  is  approximately  10.  These  two  Hj  values,  obtained  for 
each  channel,  define  the  usefiil  tangent  height  ranges  of  the  SPIRIT  HI  interferometer. 
The  results  are  shown  in  Table  3.  The  limb  radiance  spectra  required  in  these 
computations  were  generated  by  the  MODTRAN  co<le^  for  tangent  heights  Hj  below  50 
km.  and  by  the  SHARC  code*^  for  tangent  heights  above  50  km. 


Table  3.  Useful  tangent  height  ranges  for  the  six  interferometer  channels* 

Chan 

Wavelength 

range(pm) 

H-rfkm) 
for  CLIP 

HT(km)  for 

S/Nn,ax=10 

1 

17-28 

55 

120 

2 

2.6-4.9 

(NO  CLIP) 

70 

3 

5.8-8.9 

30 

90 

4(OPEN) 

4-29 

60 

160 

5 

10.4-13.2 

8 

90 

6(PREW) 

3-29 

45 

130 

*  Preliminary  estimates 


APPENDIX 

Reproduced  here  is  the  camera-ready  manuscript  of  the  paper  "Mid-Course  Space 
Experiment  (MSX):  Capabilities  of  the  LWIR  interferometer  for  remote  sensing  of  trace 
constituents  in  the  stratosphere  and  mesosphere."  The  paper  was  presented  at  the 
Atmospheric  Propagation  and  Remote  Sensing  III  conference,  part  of  SPIE's  Optical 
Engineering/Aerospace  Sensing  symposium  held  on  4-8  April  1994  at  Orlando,  FL.  The 
paper  will  appear  in  Proc.  SPIE,  vol.  2222(1994). 

The  work  reported  in  the  paper  is  ba.sed  on  the  described  noise  mcxlel  except  that 
measured  NESRs  were  used  to  scale  those  predicted  by  the  model,  as  noted  in  the  section 
Channel  NESRs. 
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ABSTRACT 

The  planned  Mid-Course  Space  Experiment  (MSX)  observations  will  include  two  experiments  for  remote  detection  of 
atmospheric  trace  constituents  above  10  km  altitude,  based  on  measurements  of  limb  spectral  radiance  by  the  cryogenic  infrared 
interferometer  and  the  ultraviolet  and  visible  spectrographic  imagers.  Species  to  be  monitored  include:  NO  and  CO-,  in  the 
thermosphere,  03and  HiO  in  the  mesosphere  and  stralo.spbere,  HNO3,  CFC-Il,  CFC-12,  N^O  and  CH4  in  the  stratosphere,  and 
CO  in  the  upper  troposphere  and  lower  stratosphere.  Quantification  of  the  altitude  profiles  of  these  species  will  give  insight  into 
processes  affecting  their  global  distributions  arid  tbt  atnursphere's  respoase  to  anthropogenic  perturbations,  and  contribute  to 
further  understanding  of  global  change.  The  timing  of  the  measurements  is  particularly  advantageous  since  they  will  likely  be 
the  only  regular  liiob  observatioas  of  trace  coastituenLs  during  the  operational  lifetime  of  the  MSX  satellite.  The  SPIRIT  III 
interfemnKter  has  a  maximum  spectral  resolution  of  1  cm  '  in  six  spectrally  isolated  channels  whose  vertical  fields  of  view  are 
between  4  and  13  km  in  line-of-sight  tangent  altitude.  The  six  channels  will  provide  spectra  over  wavelengths  in  the  2.6-28  pm 
range  for  tangent  heights  up  to  180  km.  The  capabilities  of  the  interferometer  for  the  planned  re  mote -seasing  experiments,  based 
on  predicted  instrument  noi.se  and  saturation  levels,  are  described  in  this  paper . 

I.  INTRODIKTION 

The  Mid-Course  Space  Experiment  (MSX)  is  a  pnigram  funded  and  managed  by  the  Ballistic  Missile  E)efeasc 
Orgenization.'  The  comprehensive  MSX  mi.ssion  and  suite  of  instruments  aboard  the  MSX  satellite  arc  described  by  Mill, 
ct.al.^'^  Planned  MSX  experiments^-^  will  gather  spectroscopic  and  radiometric  data  on  terrestrial,  earthlimb  and  celestial 
scenes  in  support  of  ballistic  mi.ssile  defense  objectives  as  well  as  for  .scientific  pursuits  in  atmospheric  remote  sensing  and 
astronomy. 

The  Fourier  transform  spectrometer  and  the  LWIR  radiometer  (together  called  SPIRIT  III)  collect  data  through  a  shared 
393  enr  aperture.  A  solid  hydrogen  cryogen  with  estimated  18-month  lifetime  cools  both  sensors,  their  foieoptics  and  the 
telescope  baffle.  An  additional  suite  of  instruments  (called  UVISI)  includes  five  ultraviolet  and  visible  spectrographic  imagers. 
Experiments  to  define  the  earthlimb  background  agairAt  which  midcourse  surveillance  systems  will  operate  use  the  radiometer 
as  the  primary  sensor  and  UVIS!  as  a  supporting  .vystem  to  provide  diagnostic  infomiation  on  atmospheric  composition  and 
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temperature.  They  will  measure  radiance  levels  and  rmc-scale  structure  associated  with  aurora,  airglow,  mesospheric  and 
noctiluceni  clouds,  joulc-bcated  atmospheres  and  stratospheric  warmings,  over  a  wide  range  of  viewing  conditions.  TlTie 
interferometer  will  lend  additional  diagnostic  support  by  providing  spectrally  resolved  data  over  the  broad  radiometer  bands, 
although  at  much  coarser  spatial  resolution  than  the  ladiometer. 

Two  planned  MSX  experiments  depend  on  the  capabilities  of  the  SPIRIT  III  interferometer  to  remotely  measure 
concentratioas  of  trace  species  throughout  the  stratosphere  and  above,  up  to  approximately  180  km.  Since  the  species  are 
observed  in  emission,  their  concentrations  can  be  obtained  for  a  wide  range  of  geographic  locatioas  for  both  day  and  night 
conditions.  Interferometer  measurement  programs  based  on  solar  cccultation  (like  ATMOS^)  are  more  limited  in  their  global 
coverage  and  provide  data  only  for  daytiiiK  conditions.  The  cryogenic  sensor  technology  represented  by  the  SPIRIT  III 
interferometer  has  been  demonstrated  by  the  similar  CIRRIS  1 A  interferometer  .system^  aboard  Space  Shuttle  Discovery  during 
the  eight-day  STS'39  mission. 

In  the  following  two  elections  we  de.scribc  the  interferometer's  basic  observational  capability  (its  spectral  ranges,  spectral  and 
spatial  resolutions,  and  gain  levels),  and  its  performance  (expected  noise  and  saturation  levels).  Section  4  discusses  the  two 
remote  sensing  experiments  and  some  features  of  the  automated  data  reduction  and  cataloging  procedures  that  peitain  to  ail 
experiments  using  the  interferometer.  In  Section  5  we  give  an  estimate  of  the  basic  sensing  capabilities  of  the  instrument  by 
listing  the  altitude  ranges  over  which  the  following  species  can  be  detected;  NO,  NO,,  N,0,  hN03,  CH^,  CO,  CFC-t  > ,  CFC-12, 
O^,  aerosols,  HiO  and  CO,. 

The  MSX  system  will  be  launched  in  November  1994  and  will  supply  important  data  during  a  critical  gap  of  six  or  more 
years  between  UARS'  and  the  first  EOS  environmental  satellite  (POEM-ENVISAT*^),  expected  to  ice  launched  after  the  turn  of 
the  cci  tury.  The  CLAES  and  ISAMS  sensors  on  UARS  have  ceased  operation,  while  the  currently  operating  HALOE  system 
perfoniLs  occupation  measurements  limited  to  sunrise  and  suaset,  rib  degrees  latitude,  and  cannot  measure  CFC's.  The  MLS 
sensor  on  UAR.S  can  measure  only  O3,  CIO  and  H,0,.  In  particular,  the  MSX  interferometer,  during  it.s  operational  period,  will 
be  the  only  sysferr  capable  of  measuring  vertical  profiles  of  the  critically  important  species  N,0,  HNO3,  CO,  CFC-1 1,  and  CFC- 
12  on  a  global  scale.  The  data  ba.se  of  calibrated  .MSX  nica.snremenLs  and  results  of  their  analysis  will  be  available  to  the 
scienfific  r  ommunity. 


2.  INTERFEROMETER  DESCRimON 


APPROXIMATE  distance  (km) 


-?0 


ao 


The  MSX  SPIRIT  III  interferometer  has  six  arsenic-doped  silicon  blocked  impurity  band  (Si;A.s/BIB)  detectors  with 
transimpedance  amplifiers  (TIAs)  operating  at  approximately  II  degrees  K.  Figure  1,  an  objccl-spacc  projection  of  the 
interferometer's  focal  plane,  indicates  the  sizes  and  offsets  of  the  six 
instantaneous  fields  of  view  (IFOVs).  Distances  read  from  the  right- 
hand  and  upper  scales  arc  in  a  plane  perpendicular  to  the  optic  axi.s  at 
the  line-of-sight  tangent  altitude  (LOSTA),  and  are  calculated  for 
LOSTA  =  60  km,  a  representative  value  for  the  remote  seasing 
experiments.  Optical  fillers  over  detectors  1,  2,  3,  5  and  6  limit  their 
range  of  spectral  response.  The  unfiltcred  channel  4  uses  the  full 
range  of  its  detector,  approximately  3  to  28  pni.  The  "pre whitened" 
filter  in  channel  6  passes  nearly  the  same  broad  range  of  wavelengths 
but  has  a  notch  that  attenuates  over  the  region  of  the  15-pm  CO, 
band.  For  limb  observations  the  notch  typically  effects  a  slight 
reduction  in  photon  noise  and  allows  operation  at  higher  gain,  which 
reduces  A/D  noise.  It  lowers  by  approximately  U)  km  the  LOSTA  at 
which  channel  6  is  saturated  compared  to  channel  4. 
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The  spectral  coverage  of  each  dclcctor-channcl,  its  vc:1ical  IFOV 
size  and  usable  range  in  LOSTA  are  listed  in  Table  1.  Near  the  Figure.  I.  The  interferometer  f<Kal  plane  projected  into 
bottom  of  the  usable  LOSTA  range  the  instrument  will  be  clase  to  object  space.  Distances  on  the  righi-hand  and  bottom 

saturation  (and  provide  maximum  signal  to  noise),  while  at  the  lop  scales  arc  fora  line-of-sighi  tangent  altitude  of  60  km. 

the  S/N  at  the  wavelength  of  maximum  limb  spectral  radiance  will  be 
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approximately  10.  These  tangent  height  limits  are  based  on  predictions  of  the  spectrally  resolved  earthlimb  radiance  profile  and 
current  estimates  of  the  instrument's  performance,  and  will  likely  be  refined  when  preflight  testing  and  calibration  of  the 
interferometer  are  completed.  Useful  data  at  LOSTAs  above  the  ranges  in  Tabic  1  can  be  obtained  by  coadding  interferograms, 
which  is  an  option  available  in  most  MSX  experiments. 


Tabic  1.  Characteristics  of  the  six  interferometer  channels 


Channel 

Spectral  range  (pm) 

Vertical  resolution 
(mr)  (km) 

Approximate 
LOSTA  range  (km)* 

5 

10.5-13 

1.17 

4.0 

8-90 

2 

2.6-4.9 

1.85 

6.3 

0  (or  nadir)  to  70 

3 

5.8-8.9 

3.7 

12.5 

30-90 

1 

17-28 

3.7 

12.5 

55-120 

4  (Open) 

3-28 

1.85 

6.3 

60-160 

6  (Prewhitened) 

3-26 

1.85 

6.3 

45-130 

*  Preliminary  estimates 


Sensing  of  stratospheric  trace  constituents  will  utilize  the  capabilities  of  channels  5,  2  and  3,  which  are  currently  estimated 
to  saturate  at  tangent  altitudes  of  8,  0  and  30  km,  respectively.  Channel  5  will  observe  the  emissions  of  chlorofluorocarbons, 
nitric  acid  and  other  species.  Channels  4  and  6,  with  saturation  altitudes  of  60  and  45  km,  respectively,  are  used  in  one  of  the 
two  remote  sensing  experiments,  as  described  in  Section  4.  Channel  2  will  not  saturate  even  when  pointed  BTH  (Below-The- 
Horizon)  except  possibly  when  the  BTH  scene  includes  specularly  reflected  sunlight.  Its  2.6-4.9  pm  spectral  range  provides 
support  to  SPIRIT  III  radiometer  experiments  that  will  measure  P  v-ATH  (Above-The-Horizon)  and  BTH  fine-scale  spatial 
structure  in  the  4.3  pm  CO,  region’,  while  allowing  observation  of  the  spectrum  of  sun-lit  aerosols  and  remote  sensing  of  some 
trace  species  net  observable  in  the  other  interferometer  channels  at  low  LOSTAs.  Channel  3  provides  a  capability  to  measure 
HiO,  CH4  and  NO2  concentrations  throughout  the  middle  stratosphere  and  above. 

Effects  of  nonlinearities  associated  with  the  high-performance  Si:As/TIA  combination  are  mitigated  by  recording  and 
processing  double-sided  rather  than  single-sided  interferograms.  The  interferograms  have  progranunable  lengths  equivalent  to 
unapodized  resolutions  of  roughly  1,  2  and  10  cm'*  in  wavenumber.  The  times  required  to  record  them  are  4.3,  2.2  and  0.43  s, 
respectively.  The  resolutions  of  the  device  after  minimal  apodization  to  suppress  sidelobes  are  approximately  2,  4  and  20  cm'*. 
Experiments  that  monitor  atmospheric  trace  constituents  will  produce  spectra  in  the  highest  resolution  mode. 

The  interferometer  (IFR)  system  has  a  12-bit  A/D  converter  and  seven  gain  levels  which  differ  by  multiples  of  two.  Hence,  it 
provides  an  effective  range  of  18  bits  (•  2.6  x  10^)  in  the  digitized  signal.  Users  may  program  the  gain  level  or  opt  to  use 
automatic  gain  ranging.  The  auto-gain  ranging  logic  is  described  in  the  following  section.  Additional  information  on  the 
interferometer  system  arc  given  in  the  SPIRIT  III  Seasor  User's  Guide.’** 


3.  INSTRUMENT  PERFORMANCE 

The  primary  instrument  effects  in  the  IFR  data  are  preamp- Johnson  (dark)  noise,  digitization  or  "A/D"  noise,  and  photon 
noise.  The.sc  have  corresponding  components  in  the  spectrum  obtained  by  Fourier  analysis  of  the  interferogram.  The  noise- 
equivalent  spectral  radiance  (NESR)  is  the  dark  spectrum  component,  which  dominates  when  the  scene  is  "dim",  i.e.,  at  very 
high  LOSTA.  The  lower  part  of  Fig.  2  (part  a)  shows  how  the  three  noi.se  components  vary,  at  12  pm  wavelength  in  channel  5, 
with  increasing  "band  photon  radiance"  (BPR).  BPR  is  the  photon  radiance  collected  over  the  channel's  entire  spectral  range, 
(after  accounting  for  optics  net  spectral  transnii.ssion  and  detector  quantum  efficiency),  and  determines  the  amplitude  of  the 
central  fringes  in  the  recorded  interferogram.  Photon  noise  increases  as  the  square  root  of  BPR,  while  A/D  noi.se  increa.ses 
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directly  with  each  factor-of-two  reduction  in  gain.  The  three  noise  spectral  radiances  curves  shown  in  part  (a)  are  equal  to  one- 
tenth  of  the  actual  computed  spectra  so  that  they  can  be  clearly  distinguished  from  the  curves  in  part  (b). 
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In  automatic  mode  the  gain  is  increased  (or  decreased)  when  the 
interferogram's  maximum  value  falls  below  (or  rises  above)  certain 
threshold  values.  The  new  gain  setting  applies  to  the  next 
interferogram.  Different  thresholds  for  gain  increase  and  decrease 
provide  the  hysteresis  needed  for  stability.  Thus,  the  A/D  noise  can 
have  one  of  two  levels  over  a  wide  range  of  BPR,  as  indicated  in 
Fig.  2,  part  (a).  When  the  BPR  continues  to  increase  while  the  gain 
is  at  its  lowest  value,  the  interferogram's  central  fringes  are  eventually 
clipped;  this  occurs  at  the  bottom  end  of  the  LOSTA  ranges  given  in 
Table  1.  Part  (b)  of  Fig.  2  shows  the  typical  behavior  of  spectrum 
total  noi.se,  equal  to  the  root-sum-square  of  the  three  components. 

Figure  3  shows  preliminary,  measured  NESRs  for  the  six 
interferometer  channels  for  the  highest  resolution  mode.  The  NESRs 
are  lower  for  the  other  two  resolution  modes  (and  are  reduced  when 
the  interferogram  is  apodized).  The  dashed  portions  of  the  curves  for 
channels  3,  4  and  6  are  extrapolated  from  the  data;  the  "hump"  in  the 
NESR  for  channel  6  results  from  the  prewhitening  notch.  The  other 
two  noise  components  (A/D  and  photon  noise),  being  dependent  on 
the  observed  scene,  must  be  estimated  case-by-case.  A  model  of  the 
interferometer  system  was  developed,  as  a  computer  code,  to  obtain 
the  spectrum  of  total  noise  given  the  NESR  and  a  prediction  of  the 
scene  spectral  radiance.  Since  total  noise  and  NESR  have  the  .same 
spectral  shape,  the  computation  amounts  to  a  scaling  operation.  The 
code  calculates  the  auto-gain  level  and  resultant  A/D  noise  contribution,  and  determines  whether  the  interferogram  will  be 
clipped.  It  al.so  produces  a  version  of  the  given  scene  spectrum  that  has  the  IFR  instrument  line  shape  corresponding  to  a 
specified  interferogram  length  and  apodization  function.  It  can  be  used  to  obtain  a  prediction  of  the  NESR,  ba.sed  on  data 
describing  the  system's  optical  components  and  electronics.  An  example  of  total  spectrum  noise  is  given  in  the  following 
section. 


BAND  PHOTON  RADIANCE  (photons/s  cm'  sf) 

Figure  2.  (a)  The  three  components  of  total  noise 
(times  0.1)  versus  "band  photon  radiance"  for  channel  5 
at  X.=  12  pm  and  for  1  cm  *  resolution,  (b)  The 
resultant  total  spectrum  noise  at  X  =  12  pm. 


Figure  3.  Preliminary  measurements  of  the  interferometer's  six  NESRs  for  the  highest  resolution  mode, 

unapodized  (1  cm'*  resolution  in  wavenumber). 


4.  PLANNED  EXPERIMENTS 


The  planned  MSX  Earth  Limb  Experiments  ELE-19  and  ELE-10,  entitled  "Stratospheric  Trace  Gas  Survey"  and  "Vertical 
Profile  Survey,"  respectively,  arc  designed  to  use  the  SPIRIT  III  interferometer  (and  the  UVISI  imaging  spectrographs)  to 
develop  a  "climatology"  of  earthlimb  spectral  radiance,  to  validate  physical  models  of  earthlimb  radiance,  and  to  recover  vertical 
profiles  of  radiatively  active  species.  Both  experiments  will  scan  the  limb  repetitively  in  a  sawtooth  fashion  by  moving  the  IFR's 
optic  axis  slowly  from  the  lowest  to  highest  tangent  height  and  resetting  rapidly  back  to  the  starting  height;  primary  data  is 
obtained  only  in  the  upward  scan.  Scanning  the  limb  upward  rather  than  downward  implies  lower  gain,  i.e.,  higher  A/D  noise, 
but  avoids  a  significant  loss  of  data  (repeated  clipping  of  the  interferogram's  central  fringes  prior  to  an  automatic  gain  reduction 
in  the  downward  scan).  The  upward  scan  rate  will  be  on  the  order  of  one  km/s.  TTie  resulting  4.3-km  change  in  LOSTA  as  each 
long  interferogram  is  recorded  is  found  through  simulalioas  to  produce  minimal  degradation  of  dominant  emission  features  in 
the  spectrum.  ELEi-19  will  cover  tangent  altitudes  from  6  to  70  km  and  use  primarily  data  from  IFTl  channels  5,  2  and  3.  ELE- 
10  will  scan  over  the  30-150  km  range  and  use  all  channels  of  the  IFR.  Each  data  collection  event  (DCE)  will  consist  of  a  series 
of  these  sawtooth  scans,  providing  data  over  a  wide  range  of  latitudes  and  seasons.  The  DCEs  will  be  repeated  to  generate  useful 
statistics  on  the  variability  of  trace  gas  concentrations. 

Interferometer  data  from  all  experiments  that  use  the  IFR  as  the  primary  sensor  or  in  a  support  role  will  undergo  regular, 
automated  analysis  at  the  Phillips  Laboratory  Data  Analysis  Center  as  the  data  is  received."*  The  purpose  is  to  reduce  the 
interferograms  to  calibrated  spectra,  provide  for  timely  asse.ssment  of  experiment  pcrfonnance.  catalog  the  data,  and  generate 
catalog  indices  that  flag  particular  phenomenologies  or  "events"  for  later  review  in  the  so-called  interactive  analysis.  The  events 
to  be  indexed,  which  have  easily  recognised  effects  in  the  spcctmm,  include  auroral  activity  and  stratospheric  ozone  depletion. 
Associated  with  each  arc  a  defined  set  of  calculations  that  result  in  one  or  more  "indicator"  values,  such  as  the  factor  by  which  a 
particular  integrated  band  radiance  exceeds  a  predicted  quie.sccnt  value  due  to  possible  auroral  enhancement  during  a  period  of 
high  Kp.  The  indicator  values,  corresponding  thresholds,  and  uncertainties  in  the  indicator  values  computed  from  the  noise 
model  described  earlier,  are  used  to  gauge  the  significance  of  the  event  and  assign  an  index  value  to  it.  Generally,  any  large 
differences  between  measured  and  expected  limb  radiance  spectra,  where  the  latter  are  defined  by  a  precomputed  "template" 
library,  will  be  flagged  by  a  high  index  value.  The  templates  for  the  IFR  are  maintained  as  interferograms  so  that  they  can  be 
processed  into  spectra  (as  needed,  during  the  automated  prwessing)  in  the  .same  manner  as  the  data  interferograms. 

Retrieval  of  species  vertical  profiles  from  the  reduced  spectra  imposes  a  much  greater  demand  on  computer  resources  than 
these  data  reduction  and  cataloging  procedures,  and  will  be  part  of  the  subsequent  interactive  analysis.  Some  interactive  analyses 
of  IFR  spectra  will  use  inferences  drawn  from  supporting  data  provided  by  the  UVISI  system  and  from  corollary  measurenKnfs 
obtained  by  cooperating  ground  sites  (e  g.,  ARM)  and  operational  weather  .satellites  (c.g.,  NOAA  1 1  and  12). 

Figures  4  and  5  show  typical  computer-display  color  images  (reproduced  here  in  grey-scale)  that  will  U;  produced  as 
graphical  display  products  for  ELE-19;  the  displays  for  ELE-IO  and  other  earthlimb  experiments  are  similar.  Information 
conveyed  in  these  figures  will  be  contained  in  files  created  in  the  automated  processing,  but  is  not  displayed  as  a  screen  image 
until  an  investigator  requests  the  display  in  the  interactive  review  proce.ss.  Figure  4  shows  the  measured  spectra  (actually 
template  spectra  in  this  illustration)  when  the  IFR's  optic  axis  is  at  76.8  km.  Each  of  the  six  figure  panels,  according  to  an  option 
.selected  by  the  interactive  user,  includes  the  estimated  total  noise  spectrum  and  shows  the  measured  limb  spectral  radiance  only 
where  it  exceeds  the  total  noise.  Figures  4  and  5  represent  concepts;  codes  to  produce  them  do  not  yet  exist. 

The  top-central  portion  of  the  display  in  Fig.  4  gives  basic  date,  time  and  pointing  information,  and  shows  options  .selected 
in  the  principal  investig.ator's  expi'rimeni  script,  such  as  the  apodi/ation  function,  which  affects  the  spectral  resolution  and  levels 
of  sidelobes  and  noise  in  the  displayed  spectra.  The  IFOV  map  at  the  upper  left  is  ccxled  in  six  colors  to  correspond  to  the  six 
displayed  spectra.  One  can  riKsuse-click  on  "buttons"  in  the  upper  right  to  view  a  new  display  showing  the  spectrum  for  only  one 
channel  on  an  expanded  scale,  and  to  sujxrimposc  the  appropriate  zodiacal  background  (ZOD)  spectrum,  an  estimate  of  Non- 
Rejected  Earth  Spectral  Radiance  (NRESR)  or  the  template  spectrum;  the  latter  three  are  generated  routinely  in  the  automated 
processing  as  each  data  interferogram  is  transformed  to  a  spectrum.  Figure  5  shows  the  result  of  a  user's  decision  to  display  only 
the  channcl-3  spectrum,  with  ZOD,  NRESR  and  template  sjx-ctra  supenniposcd.  The  upper  panel  indicates  the  location  of  the 
displayed  spectrum  and  others  measured  in  the  ELE-19  sawtooth  limb  scan.  Other  interactive  options  in  the  two  displays  allow 
stepping  over  sei|Ueiitial  spectrum  sets  within  the  current  DCE  or  advancing  to  the  next  DCE. 
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5.  PRKDK'TEI)  SENSING  CAPABILITIES 


Detailed  computer  simulations  will  be  performed  to  quantify  the  expected  precision  and  accuracy  in  species  vertical  profiles 
to  be  recovered  front  ELE-19  and  ELE-10  data.  Until  these  are  completed,  our  estimates  of  remote  sensing  capability  are  based 
on  calculatioas  of  "species  contributions".  In  this  context  a  contribution  is  defined  as  the  effect  a  particular  species  has  in  the 
observed  spectrum  when  introduced  to  an  atmosphere  that  contains  all  but  this  species.  By  intercomparing  the  relative 
contributions  of  all  expected  species  and  total  noise,  for  all  IFR  channels  over  a  range  of  LOST As,  one  can  at  least  e.stimatc 
vertical  ranges  for  the  recovered  profiles,  if  not  accuracies  in  the  species  concentrations.  This  procedure,  admittedly,  is  not  as 
reliable  or  informative  as  the  planned  full-blown  simulation  and  error  analysis  in  which  retrieval  algorithms,  like  those  described 
by  Marks  and  Rodgers'*,  will  be  applied  to  noise-contaminated  .synthetic  data. 

Figure  6  shows  the  total  spectral  radiance,  predicted  total  system  noise  level  and  all  significant  .species  contributions  over  the 
wavelength  range  of  channel  5  when  its  detector  is  centered  on  LOSTA  =  15  km.  The  noise  level  is  based  on  early  estimates  of 
pcrfonnance  and  is  expected  to  be  somewhat  lower  when  recalculated  from  final  test  data.  The  spectral  resolution,  corresponding 
to  triangular  apodization,  is  a  constant  2  cm"'  in  wavenumber.  Calculation  of  these  species  contributions  and  total  radiance  used 
a  version  of  the  MODTRAJM  Code'*  modified  by  the  Philii|>s  Laboratory  to  include  CCLFt  (CFC-12  or  F-12),  CCI3F  (CFC-1 1  or 
F  ll)  and  CCI4.  It  is  seen  that  these  compounds,  as  well  as  nitric  acid,  ozone  and  COt,  have  contributions  at  LOSTA  =15  km 
that  do  not  seriously  mask  each  other  and  arc  well  above  both  the  aerosol  thermal  background  and  sensor  total  noi.se  level.  The 
procedure  to  be  used  for  the  limb  retrievals  will  attempt  to  recover  these  species  simultaneously,  using  all  wavelengths  in  the  data 
spectrum. 


Figure  6.  Species  contributions  and  total  noise  computed  for  channel  5  fora  tangent  height  of  15  km.  The  spectral  signatures 
and  interferometer  noi.se  level  are  early  estimates  corresponding  to  a  spectral  resolution  of  2  cm  '. 

it  was  determined  from  similar  calculations  that  the  interferometer  will  be  able  to  remotely  sense  the  species/parameters 
appearing  in  Fig.  7,  ove."^  the  altitude  ranges  shown.  Note  that  many  of  the  species  arc  detectable  down  to  approximately  10  km, 
the  LOSTA  at  which  channel  5  is  expected  to  be  nearly  saturated.  The  ultraviolet  and  visible  instruments  enhance  the  indicated 
IFR  capability  by  adding  some  species  and  extending  (during  daytime)  the  ranges  shown  in  Fig  7  for  NO,  and  aerosols.'^ 
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Figure  7,  Measurement  capabilities  of  the  MSX  SPIRIT  III  interferometer.  The  ordinate  scale  gives 
preliminary  estimates  of  the  tangent  heights  over  which  species  vertical  profiles  can  be  retrieved. 


6.  SUMMARY 

The  comprehensive  mission  of  the  MSX  program  and  its  extensive  instrumentation  suits  it  well  to  the  study  of  scientific 
issues  having  global  significance,  such  as  changes  in  the  distributions  of  atmospheric  trace  species,  'fhe  planned  earthlimb 
observations  include  two  interferometer  experiments  dedicated  to  measuring  vertical  profiles  of  NO,  NO2  ,  NjO,  HNO3,  CH^, 
CO,  CFC-11,  CFC-12,  O3,  aerosols,  H2O  and  CO2  over  the  altitude  ranges  shown  in  Fig.  7.  These  results  will  be  supported  by 
inferences  drawn  from  MSX  ultraviolet  and  visible  measurements,  and  by  data  from  corollary  ground  sites  and  weather  satellites. 
The  set  of  measured  spectra,  to  be  acquired  over  an  18-montb  period  when  no  comparable  data  will  be  available  from  other 
satellite  sensors,  will  give  insight  into  chemical  and  dynamical  processes  affecting  the  global  distributions  of  the  above  species, 
and  nrovide  evidence  of  short  term  changes  in  chemistry  and  climate. 
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